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INTRODUCTION 


The  existence  of  trace  amounts  of  atomic  sodium  (Na)  In  the  Earth's  upper 
atmosphere  was  deduced  In  the  earlier  part  of  this  century  from  observations 
of  nighttime  airglow  in  the  D  lines  of  Na  (Chapman,  1939,  and  references  there 
given).  Subsequent  measurements  of  resonance  radiation  in  the  D  lines  were 
carried  out  at  twilight  and  during  the  day  by  numerous  investigators  in  the 
period  prior  to  1970  (see  review  by  Kvifte,  1973).  These  measurements  deduced 
the  concentration  of  atomic  sodium  in  a  narrow  layer  centered  at  around  90  km. 
The  advent  of  laser  radar  techniques  (lidar)  has  produced  a  wealth  of  infor¬ 
mation  on  the  altitude  distribution  and  the  diurnal  and  seasonal  variation  of 
the  sodium  layer  during  the  past  15  years  (Gibson  and  Sandford,  1971,  1972; 
Megie  and  Blamont,  1977;  Simonich  et  al.,  1979;  Clemesha  et  al.,  1982;  Granier 
and  Megie,  1982).  Although  different  mechanisms  were  originally  proposed  to 
explain  the  origin  of  this  sodium  layer  (Kvifte,  1973),  experimental  evidence 
and  present  theoretical  models  favor  a  source  from  ablation  of  meteorites  oc- 
curing  between  80  and  100  km  (Megie  and  Blamont,  1977;  Gadsden,  1968;  Hunten 
et  al.,  1980). 

The  first  photochemical  model  for  sodium  compounds  was  proposed  in  1939 
by  Chapman  to  explain  the  nighttime  D-line  emission.  This  mechanism  proceeds 
through  the  oxidation-reduction  cycle  involving  Na  and  sodium  monoxide  (NaO) 

Na  +  O3  -*■  NaO  +  0  (RR1)* 

NaO  +  0  +  Na(2P,2S)  +  02 ,  (RR2) 

where  the  dissociation  energy  of  02  is  partially  channelled  into  the  excita¬ 
tion  of  Na  to  the  2P  state,  with  subsequent  radiation  in  the  D  lines.  The 
validity  of  the  above  mechanism  has  been  confirmed  by  kinetic  studies  and  at¬ 
mospheric  observations  during  the  past  40  years.  The  evidence  against  the 
possibility  of  other  mechanisms  has  been  recently  summarized  by  Bates  and  Ojha 
(1980). 

The  existence  and  chemistry  of  sodium  compounds  other  than  atomic  sodium 
is  more  uncertain,  because  of  a  lack  of  laboratory  kinetic  data  and  atmospheric 


*RR  numbers  refer  to  reaction  numbers  in  Table  1 


measurements.  The  non  proton-hydrate  positive  ions  observed  in  the  middle 
atmosphere  by  Arnold  et  al.  (1978)  were  originally  Interpreted  to  be  of  the 
form  Na0H2+(Na0H)m(H20)n  (Ferguson,  1978).  This  interpretation  led  Liu  and 
Reid  (1979)  to  postulate  the  production  of  NaOH  through  reactions  such  as 


Na02  +  OH  NaOH  +  02, 


Na02  +  H  NaOH  +  0. 


where  Na02  is  produced  by  the  reaction 


Na  +  02  +  M  ♦  Na02  +  M. 


(RR10) 


(RR1 1 ) 


(RR3) 


This  model  predicted  rapid  conversion  of  Na  to  NaOH  below  an  altitude  of 
90  km.  Direct  production  of  NaOH  from  NaO  was  proposed  by  Murad  and  Swider 
(1979)  through  the  reaction 


NaO  +  H20  -*•  NaOH  +  OH. 


(RR7) 


Although  more  recent  measurements  of  stratospheric  positive  ions  have  elimi¬ 
nated  the  possibility  of  a  sodium  core  (Arnold  et  al.,  1981),  the  NaOH  sources 
represented  by  reactions  (RR10.RR11)  and  (RR7)  still  remain  a  possibility. 

The  abundance  of  NaOH  is  also  of  interest  for  stratospheric  ozone  chemis¬ 
try.  Murad  and  Swider  (1979)  and  Murad,  Swider,  and  Benson  (1981)  pointed  out 
that  the  reaction 

NaOH  +  HC1  +  NaCl  +  H2O  (1) 

could  constitute  an  Important  reaction  pathway  for  stratospheric  HC1  if  its 
rate  were  of  the  order  of  10~12  cm2  s”*,  and  if  the  NaOH  concentration  were 
~  105  cm-3.  Furthermore,  if  NaCl  were  sufficiently  stable  to  act  as  a  nucle¬ 
ating  agent  and  be  incorporated  into  aerosols,  reaction  (1)  could  provide  an 
important  sink  for  stratospheric  chlorine.  The  value  of  the  rate  for  reaction 
(l)  has  recently  been  measured  to  be  2.5  x  10”1®  cm3  s-*  (Silver,  Stanton, 
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Zahniser,  and  Kolb,  1984a).  However,  Rowland  and  Rogers  (1982)  indicated  that 
photolysis  of  NaCl 

NaCl  +  hv  Na  +C1  (2) 

could  have  a  value  of  ~  10~2  s'"1  if  the  observed  UV  absorption  cross  section 
of  NaCl  is  taken  to  be  its  photodissociation  cross  section.  In  this  case,  the 
sequence  of  reactions  (1)  and  (2)  acts  as  a  catalytic  cycle  converting  the 
relatively  stable  HC1  into  active  Cl.  Determination  of  the  importance  of  re¬ 
actions  (1)  and  (2)  in  the  stratosphere  awaits  the  measurement  of  unknown 
kinetic  rates,  in  particular  the  production  of  NaOH  (RR7). 

More  recently,  the  rate  for  RR3  has  been  measured  (Husain  and  Plane, 

1982;  Silver,  Zahniser,  Stanton,  and  Kolb,  1984b)  to  have  a  value  three  orders 
of  magnitude  higher  than  that  used  in  Liu  and  Reid's  model.  Such  high  values 
for  the  rate  of  RR3  raised  the  possibility  of  the  existence  of  large  amounts 
of  Na02  in  the  upper  atmosphere,  with  consequently  low  concentrations  of 
atomic  sodium.  Analysis  of  the  latest  thermocheraical  data  led  Sze,  Ko,  Murad, 
and  Swider  (1982)  to  propose  the  reaction 

Na02  +  0  NaO  +  02  (RR6) 

as  a  potentially-f ast  recycling  mechanism  for  Na02» 

Our  present  understanding  of  the  chemistry  of  sodium  compounds  is  summar¬ 
ized  in  Figure  1,  based  on  the  scheme  of  Sze,  Ko,  Murad,  and  Swider  (1982). 
Because  of  the  difficulty  in  manipulating  gas-phase  Na  in  the  laboratory, 
kinetic  data  for  reactions  involving  sodium  are  almost  nonexistant.  Besides 
the  measurement  of  the  rate  for  RR3,  experimental  data  exists  only  for  the 
photolysis  of  NaOH  (RR4)  (Rowland  and  Makide,  1982).  However,  constraints  can 
be  placed  on  the  rate  for  RRl  and  RR2  from  the  observed  D-line  emission,  the 
measured  background  concentrations  of  0  and  03,  and  kinetic  considerations 
(Bates  and  Ojha,  1980).  Additional  constraints  on  other  reaction  rates  can  be 
obtained  by  comparing  calculations  with  observations  of  the  distribution  (Liu 
and  Reid,  1979;  Sze,  Ko,  Murad,  and  Swider,  1982;  Thomas,  Isherwood,  and  Bow¬ 
man,  1983)  and  diurnal  variation  of  sodium  (Kirchoff  and  Clemesha,  1983; 
Kirchoff,  1983).  Estimated  rates  are  given  in  Table  1  where  laboratory 
surements  do  not  exist. 
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In  Figure  1,  the  Na-NaO  on  the  left-hand  side  represents  the  original 
Chapman  scheme.  The  rates  for  RR1  and  RR2  are  sufficiently  fast  that  the 
[Na]/[NaO]  should  be  proportional  to  the  [0] / [O3 ]  at  any  time  of  the  day. 

Thus,  in  the  absence  of  other  sodium  species,  the  diurnal  behavior  of  the  Ha 
and  NaO  would  be  determined  entirely  by  the  diurnal  variations  of  0  and  O3. 

The  species  Na02  and  NaOH  act  as  temporary  reservoirs  for  the  sodium  species 
during  the  diurnal  cycle.  With  the  short  photolysis  lifetime  for  NaOH,  diur¬ 
nal  variation  of  Na  is  expected  if  the  density  of  NaOH  is  sufficiently  high. 

The  effect  of  Na02  on  the  diurnal  variation  of  Na  would  depend  on  the  magni¬ 
tude  of  the  recycling  reactions  (RR5  and  RR6). 

We  present  in  this  report  the  results  of  our  calculations  of  the  diurnal 
variation  of  sodium  compounds.  Comparison  of  our  results  with  the  observed 
properties  that  are  determined  by  photochemical  processes  can  place  additional 
constraints  on  some  of  the  reaction  rates.  In  particular,  we  will  see  that 
although  RR6  must  be  slow,  a  fast  rate  for  RR7  is  consistent  with  present  ob¬ 
servations.  Our  study  will  also  isolate  the  chemical  cycles  controlling  the 
distribution  and  diurnal  behavior  of  the  sodium  layer,  thus  estimating  the 
sensitivity  of  our  results  to  unknown  kinetic  rates  and  properties  of  the 
background  atmosphere. 

Section  2  summarizes  the  experimental  information  about  atmospheric 
sodium  available  at  present.  Our  model  is  described  in  Section  3,  and  the 
results  of  our  calculations  are  presented  in  Section  4.  Conclusions  of  our 
study  are  discussed  in  Section  5. 

2.  SUMMARY  OF  EXPERIMENTAL  DATA 

An  extensive  data  base  now  exists  providing  detailed  information  on  the 
diurnal  and  seasonal  behavior  of  the  sodium  layer.  Some  of  these  properties 
are  directly  attributable  to  transport  phenomena  on  a  seasonal  (Gibson  and 
Sandford,  1971;  Megie  et  al.,  1978;  Kirchoff  et  al.,  1981;  Kirchoff  and  Cle- 
mesha,  1983)  or  diurnal  time  scale  (Megie  and  Blamont,  1977;  Clemesha  et  al., 
1978;  Kirchoff  et  al.,  1981;  Clemesha  et  al.,  1982).  Nevertheless,  we  can 
deduce  from  the  data  certain  characteristics  that  have  direct  relevance  to 
photochemical  processes.  The  experimental  data  to  be  used  in  constraining  our 
model  are  as  follows: 


a)  Total  Abundance  of  Atomic  Sodium:  Lidar  observations  have  con¬ 
strained  the  diurnal  variation  of  the  Na  column  density  to  less  than 
152  (Gibson  and  Sandford,  1972;  Clemesha  et  al.,  1982).  The  abso¬ 
lute  magnitude  of  the  total  abundance  ranges  from  a  summer  minimum 

9  —2  Q  _  n 

of  ~  2  x  10  cm  to  a  winter  maximum  of  6-8  x  107  cm  ,  yielding  an 

9  —2 

average  column  density  of  4  x  10  cm  . 

b)  Altitude  Profiles  of  Atomic  Sodium;  Lidar  observations  of  the  sodi¬ 
um  layer  at  different  latitudes  have  consistently  deduced  a  peak 
density  of  3-5  x  10  cm  at  90-92  km  during  periods  of  low  meteor- 
itic  activity  (Megie  and  Blamont,  1977;  Simonich  et  al.,  1979).  The 
peak  height  and  density  correlates  with  the  seasonal  change  in  total 
abundance,  with  high  column  densities  corresponding  to  low  peak  al¬ 
titudes  and  high  peak  densities  (Gibson  and  Sandford,  1971). 

Diurnal  variations  in  the  peak  height  have  also  been  observed,  but 
they  are  probably  due  to  transport  effects  of  gravity  waves  (Cle¬ 
mesha  et  al.,  1982). 

c)  Diurnal  Variation  of  Local  Densities  of  Atomic  Sodium:  The  densi¬ 


ties  measured  by  lidar  techniques  exhibit  a  diurnal  variation  below 
84  km  (Simonich  et  al.,  1979;  Clemesha  et  al.,  1982).  The  nocturnal 
decrease  in  sodium  densities  observed  by  these  authors  is  about  a 
factor  of  2  at  83  km  and  10  at  81  km. 

d)  D-line  Emission  Rates:  Measurements  of  the  nighttime  D-line  emis¬ 
sion  in  the  southern  hemisphere  (low  latitude)  yield  minimum  values 
of  30  Rayleigh  (R)  during  winter  and  60  to  100  R  at  equinox  (Kir- 
choff  et  al.,  1981).  Mean  intensities  in  the  northern  hemisphere 
range  from  35  to  100  R  (Fukuyama,  1977).  The  seasonal  variation  in 
the  northern  hemisphere  measurements  is  most  pronounced  at  low  lati¬ 
tudes,  with  amplitudes  in  agreement  with  the  southern  hemisphere 
results.  The  emission  rate  also  exhibits  a  small  nocturnal  varia¬ 
tion  of  30%  symmetric  around  a  minimum  at  midnight  (Kirchoff  et  al., 
1981). 

A  summary  of  the  above  experimental  constraints  is  given  i.n  Table  2.  We 
will  concentrate  on  the  diurnal  variation  of  the  above  quantities.  However, 
the  seasonal  behavior  will  be  useful  for  placing  bounds  on  the  experimental 
values  used  in  comparison  to  our  one-dimensional  model. 


3. 


DESCRIPTION  OF  MODEL 


The  one-dimensional  model  of  mesospheric  trace  gases  covers  the  altitude 
range  from  60—100  km.  It  is  an  extension  of  the  model  previously  used  in 
studies  of  stratospheric  chemistry  (Sze,  Ko,  Specht,  and  Livshits,  1980). 
Concentrations  of  short-lived  photochemical  species  are  calculated  using  the 
grouping  technique  previously  employed  in  stratospheric  models  (Sze  et  al., 
1980). 

The  species  calculated  in  our  model  include: 
oxygen  species,  0X:  O^D),  0,  O3 

odd  hydrogen  species,  H0X:  H,  OH,  HO2 ,  H2O2 

odd  nitrogen  species,  N0X:  N,  NO,  NO2,  HNO3,  NO3 ,  N2O5,  HO2NO2 
sodium  species,  NaX:  Na,  NaO,  Na02 ,  NaOH 
In  addition,  02(*A  )  is  also  included.  Although  the  NO  species  and  02(*A  ) 

1  o  X  g 

are  of  minor  importance  in  the  neutral  chemistry  of  the  odd-oxygen  and  odd- 
hydrogen  species,  they  are  included  in  the  chemistry  scheme  in  anticipation  of 
their  roles  in  the  ion  chemistry.  The  effect  of  the  odd-chlorine  and  methyl 
families  on  the  above  species  is  negligible  at  these  altitudes;  thus,  they 
have  not  been  included  in  the  model. 


Background  species  which  are  held  fixed  in  the  model  include  H2 ,  H2O,  CO, 
N2O,  O2 »  and  n2*  Densities  of  the  major  constituents  (N2  and  O2)  and  tempera¬ 
tures  are  taken  from  the  U.S.  Standard  Atmosphere  (1976).  The  densities  of 
the  background  trace  species  CO,  H2,  and  H2O,  are  taken  from  Allen,  Lunine, 
and  Yung  (1984). 

Photochemical  lifetimes  for  H0X  and  Ox  are  shorter  than  one  day  at  alti¬ 
tudes  below  80  km.  The  concentrations  of  all  constituents  are  calculated  from 
the  time-dependent  equation 


dn^ 

dT 


Pi  "  LiV 


(3) 


where  n^,  ,  P^ ,  and  L^  denote  the  density,  production,  and  loss  frequency  for 
the  ic^  species.  Equation  (3)  is  solved  subject  to  periodic  boundary  condi¬ 
tions  over  a  24-hour  period, 


n  (t  +  24  hrs.)  »  n  (t). 


(4) 


Note  that  the  effect  of  transport  is  ignored  in  this  region. 

Time  constants  for  H0X,  0X,  and  N0X  become  much  longer  than  one  day  above 
90  km.  The  densities  of  the  above  constituents  are  then  determined  by  down- 
gradient  diffusion.  Although  it  is  within  the  capability  of  the  model  to  cal¬ 
culate  the  distribution  of  the  long-lived  species,  we  have  fixed  their  densi¬ 
ties  in  our  calculations  to  the  profiles  given  by  Allen  et  al.,  (1984)  for  H0X 
and  0  ,  and  by  Solomon,  Crutzen,  and  Roble  (1982)  for  NO  . 

X  X 

Sharp  vertical  gradients  in  the  densities  of  0X  and  H0X  result  in  very 
short  (~  1  day)  transport  time  constants  between  80  and  85  km.  The  diurnal 
variation  of  the  above  species  is  then  determined  by  both  photochemical  and 
transport  processes.  The  effects  of  transport  are  approximately  accounted  for 

ic 

in  the  present  calculation  by  introducing  an  extra  time-independent  term, 
in  equation  (3).  The  magnitude  and  sign  of  P^*  is  chosen  such  that  the  solu¬ 
tions  of  equation  (3)  at  noon  agree  with  the  values  calculated  by  Allen  et 
al.,  (1984)  at  the  same  time. 

The  total  density  of  sodium  species  (NaX)  is  calculated  assuming  an  in¬ 
coming  flux  of  2  x  104  cm  2  s  1  at  100  km.  The  sodium  is  transported  by  dif¬ 
fusion  into  the  troposphere  where  it  is  removed  by  heterogeneous  processes 
with  an  effective  lifetime  of  10  days.  The  values  of  the  diffusion  coeffi¬ 
cients  below  70  km  are  the  same  as  those  used  in  the  stratospheric  model  (Sze, 
Ko,  Specht,  and  Livshits,  1980).  Above  70  km,  a  value  of  106  cm2  s_1  is 
adopted  as  suggested  by  Hunten  (1975).  The  calculated  total  NaX  volume  mixing 
ratio  ranges  from  2.3  x  10  9  at  100  km  to  1  x  10  11  at  60  km. 

The  chemical  scheme  for  the  oxygen,  hydrogen,  and  nitrogen  species  is 
similar  to  the  one  used  in  the  stratospheric  model.  Reaction  rates  are  taken 
from  the  latest  NASA  review  (JPL,  1982),  with  the  corrections  suggested  by 
Allen  et  al.,  (1984).  The  reactions  involving  sodium  species  are  taken  from 
Sze,  Ko,  Murad,  and  Swider  (1982)  as  given  in  Table  1. 

4 .  RESULTS 

Diurnal  profiles  for  the  main  0X  species  (0,  O3)  at  82,  84,  86,  and  88  km 
are  shown  in  Figures  2a,  2b,  2c,  and  2d,  respectively;  similar  profiles  at  the 
same  altitudes  are  shown  for  the  H0X  species,  H,  OH,  and  HO2  in  Figures  3a, 

3b,  3c,  and  3d.  The  calculations  are  performed  for  38°N,  summer  solstice  con¬ 
ditions.  Our  results  are  similar  to  the  diurnal  variation  calculated  by  Allen 


et  al.,  (1984).  The  Ox  and  HOx  families  are  dominated  by  0  and  H,  respective¬ 
ly.  The  change  in  the  diurnal  behavior  of  these  species  with  altitude  illus¬ 
trate  the  transition  between  the  photochemically-cont rolled  region  below  82  km 
and  the  diffusion-controlled  region  above  90  km.  Since  the  sodium  chemistry 
above  60  km  has  very  little  effect  on  the  background  neutral  chemistry,  the 
above  profiles  are  used  for  all  the  different  models  of  sodium  considered  be¬ 
low. 

The  results  shown  in  Figure  2  indicate  that  the  [0 ] / [O3 )  ratio  is  102  - 
103  between  80  and  90  km.  In  the  absence  of  temporary  reservoirs,  [Na]  and 
[NaO]  would  be  controlled  by  the  Chapman  scheme.  The  rates  for  RR1  and  RR2 
given  in  Table  1  would  then  imply  a  [Na]/[NaO]  ratio  of  101  -  102,  placing  the 
atomic  sodium  peak  well  below  80  km.  Thus,  the  observed  position  of  the  Na 
peak  and  the  constraints  on  RR1  and  RR2  as  suggested  by  thermal  chemical  data 
strongly  suggest  the  existence  of  temporary  reservoirs  for  atomic  sodium. 

Our  choice  of  models  for  the  sodium  chemistry  focuses  on  the  interaction 
of  the  Na-NaO  species  with  the  temporary  reservoirs,  Na02  and  NaOH.  In  the 
chemical  scheme  illustrated  in  Figure  l,  the  interaction  between  Na  and  Na02 
is  controlled  by  the  rapid  formation  of  Na02 ,  (reaction  RR3) ,  and  recycling  of 
Na02  to  Na-NaO  through  RR6.  Similarly,  coupling  to  NaOH  is  controlled  by  RR7 
and  the  photolysis  reaction  RR4.  Since  experimental  data  exists  on  both  RR3 
(Silver  et  al.,  1984b)  and  RR4  (Rowland  and  Makide,  1982),  our  study  will  con¬ 
centrate  on  reactions  RR6  and  RR7. 

Four  models  were  considered,  representing  different  strengths  of  coupling 
between  the  Na-NaO  species  and  the  Na02  and  NaOH  temporary  reservoirs.  The 
adopted  parameters  for  these  models  are  listed  in  Table  3.  Model  A  represents 
the  conditions  for  weakest  coupling  to  temporary  reservoirs  consistent  with 
measured  kinetic  rates.  The  value  adopted  for  reaction  RR6  is  small,  but  suf¬ 
ficient  to  recycle  enough  Na02  back  to  Na-NaO  so  that  Na  does  not  disappear 
cot..  ely  (Sze  et  al.,  1982)  due  to  the  fast  rates  measured  for  reaction  RR3. 

Produ  on  of  NaOH  by  RR7  is  assumed  to  be  zero.  The  coupling  between  NaO-Na 

and  N?'  Is  weak,  since  production  of  NaOH  occurs  only  via  reactions  RR10  and 
RR11  hrough  Na02  as  an  intermediate.  Model  C  represents  conditions  of  maxi¬ 
mum  cc.  tpling  to  both  temporary  reservoirs  with  rates  for  reactions  RR6  and  RR7 
near  tie  kinetic  limit.  Models  B  and  D  exhibit  combinations  of  strong  and 
weak  coupling  to  the  different  temporary  reservoirs. 


Density  (cm3)  Density  (cm 
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Table  2:  Experimental  Data  on  Sodium 


Sodium  Layer 

Peak  Altitude 
Peak  Density 

Total  Abundance 


Diurnal  Variation 


Total  Abundance3 


l Na 1  Dusk/ [Nal Dawn 


D-line  Emission  Rate 


References 


90  -  92  km 


3-5  x  10^  cm“^ 


4  x  109  cm~^ 


(2-8  x  109  cm-^) 


<  15% 


2  (83  km) 
10  (81  km) 


~  50  R 
(30-100  R) 


Megie  and  Blamont,  1977; 
Simonich  et  al.,  1979 


Gibson  and  Sandford,  1972; 
Clemesha  et  al.,  1982 


Gibson  and  Sandford,  1971; 
Clemesha  et  al.,  1982 


Simonich  et  al.,  1979; 
Clemesha  et  al.,  1982 


Fukuyama,  1977;  Kirchoff  et 
al.,  1981 


Percent  deviation  from  median  abundance 


y 


v'1 


Table  3:  Adopted  Models 


k  (Na02  +  0)  k  (NaO  +  H20) 

(cm3  s~l)  (cm3  s'1) 


Model  A 

*— 

o 

1 

u> 

0.0 

Model  B 

10"13 

io~io 

Model  C 

10_l° 

10_10 

Model  D 

10-10 

0.0 

Reaction  Scheme  and  Rate  Constants 


Reaction 

Rate  Coefficient3 

1. 

Na  +  03 

NaO  +  02 

3(-10) 

[b] 

2. 

NaO  +  0 

Na  (2P,  2S)  +  02 

4(-ll) 

tc) 

3. 

Na  +  02  +  N2 

Na02  +  N2 

1.9  X  10_30(T/300)_1 

[d] 

4. 

NaOH  +  hv 

Na  +  OH 

2(-3) 

[e] 

5. 

Na02  +  hv 

Na  +  02 

l(-4) 

[c] 

6. 

Na02  +  0 

NaO  +  O2 

K-10),  1  (—13 ) 

[c] 

7. 

NaO  +  H2O 

+ 

NaOH  +  OH 

0,  1(-10) 

[c] 

8. 

NaO  +  O3 

+ 

Na02  +  O2 

5(— 1 1 ) 

[f] 

9. 

NaO  +  H 

Na  +  OH 

1 (—14) 

[g] 

10. 

Na02  +  OH 

NaOH  +  O2 

K-ll) 

[f] 

11. 

Na02  +  H 

♦ 

NaOH  +  0 

1(-13) 

fg] 

12. 

NaOH  +  H 

*► 

Na  +  H2O 

1 -4(— 12) 

If] 

13. 

NaOH  +  (lD) 

NaO  +  OH 

K-10) 

[g] 

14. 

Na  +  H02 

♦ 

NaOH  +  0 

i(-io) 

If] 

15. 

NaO  +  O3 

Na  +  202 

K-10) 

Ig] 

16. 

NaO  +  HO2 

NaOH  +  02 

K-ll) 

[f] 

aRead  3(-10)  as  3  x  10  Units  are  cm3  s  1 ,  except  cm6  s"1  for  (3)  and 

s  1  for  (4),  (5),  and  (15).  [b],  Kolb  and  Elgin  (1976).  [c],  Sze  et  al 

(1982).  [d].  Silver  et  al.  (1984).  [e],  Rowland  and  Makide  (1982). 

[f],  Liu  and  Reid  (1979).  [g],  Kirchoff  and  Clemesha  (1983). 


change  with  02+,  with  an  effective  J  rate  of  ~  10~5  s”1  (Swider,  1970;  1984). 
Loss  of  Na+  occurs  primarily  through  clustering  with  N2» 

Na+  +  N2  +  M  *Na+  •  N2  +  M  (7) 

followed  by  dissociative  recombination  of  Na+  •  N2  (Richter  and  Sechrist, 
1979).  The  partitioning  between  Na  and  Na+  is  determined  by  the  rate  of  the 
clustering  reaction  (7),  which  has  not  been  measured.  If  one  uses  the  esti¬ 
mated  rate  (Richter  and  Sechrist,  1979)  based  on  analogous  reactions  involving 
NO*,  as  much  as  half  of  the  atomic  sodium  could  be  in  ionized  form  at  100  km. 
This  would  reduce  the  amount  of  atomic  Na  calculated  in  the  model  and  intro¬ 
duce  diurnal  variation  of  Na  at  these  altitudes  bringing  the  result  into  much 
better  agreement  with  the  measurements  shown  in  Figure  4.  Ionization  of  Na 
could  thus  account  for  the  small  Na  scale  height  (2  to  3  km)  observed  above 
the  peak  (Simonich  et  al.,  1979;  Granier  and  Megie,  1982). 

In  conclusion,  the  results  of  our  study  point  out  the  importance  of  the 
Na02  temporary  reservoir  in  determining  the  diurnal  variation  of  atomic  sodi¬ 
um.  The  diurnal  variation  of  atomic  oxygen  thus  plays  a  dominant  role  in  de¬ 
termining  diurnal  changes  in  Na.  Present  experimental  data  constrains  the 
recycling  of  Na02  through  reaction  with  0  (RR6),  limiting  the  rate  to  consi¬ 
derably  less  than  10_l°  cm3  s_1,  but  at  least  10”13  cm3  s_1 .  The  data,  how¬ 
ever,  does  not  rule  out  production  of  large  amounts  of  NaOH  through  the  reac¬ 
tion  of  NaO  with  H2O  (reaction  RR7).  Further  clarification  of  this  problem 
requires  measurements  of  unknown  reaction  rates  (in  particular,  RR4,  RR5,  RR6, 
and  RR7),  continuing  high-resolution  observations  of  sodium  species  in  the 
upper  atmosphere,  and  further  modelling  incorporating  tidal  and  gravity  wave 
effects  on  both  sodium  species  and  background  trace  gases. 


Figure  2  and  that  of  Na  in  Figure  6  suggests  general  agreement  between  both 
diurnal  changes.  In  fact,  the  ratio  t^Dawn  is  ^*3  at  82  km  and  1.7 

at  84  km,  similar  to  the  sodium  ratio  in  Model  A.  The  additional  coupling  to 
NaOH  is  reflected  in  larger  ratios  for  Models  B  and  C.  Since  the  sodium  peak 
is  at  82  km  for  Model  D,  NaC>2  is  no  longer  the  major  sodium  species,  and  the 
local  nocturnal  variation  is  smaller  than  for  the  previous  models. 

Given  that  diurnal  variation  of  atomic  sodium  depends  crucially  on  the 
diurnal  variation  of  atomic  oxygen,  the  observed  absence  of  appreciable  diur¬ 
nal  change  in  the  total  Na  abundance  requires  that  the  bulk  of  the  sodium  be 
located  in  a  region  where  0  is  constant  throughout  the  day,  i.e.,  about  85 
km.  Once  the  total  NaX  is  fixed,  the  position  of  the  peak  is  determined  to 
first  order  by  the  magnitude  of  RR3.  Lowering  the  peak  by  modifying  the  NaX 
profile  through  changes  in  the  Na  source  or  eddy  diffusion  coefficient  could 
be  reflected  in  a  larger  diurnal  variation  of  the  column  density  of  sodium. 
Since  the  time  constant  for  the  Na-Na02  exchange  ranges  from  seconds  to  a  few 
minutes  below  90  km,  any  modelling  of  transport  effects  (i.e.,  tides,  gravity 
waves)  must  include  a  self-consistent  treatment  of  the  atomic  oxygen  variation. 

The  photolytic  reactions, 

NaOH  +  hv  -*■  Na  +  OH,  (RR4) 

Na02  +  hv  -►  Na  +  02 ,  (RR5) 

could  also  play  an  important  role  in  the  diurnal  variation  of  sodium  if  their 
values  were  very  different  from  those  assumed  here.  Reaction  RR5  could  com¬ 
pete  with  the  reaction  with  0  (RR6)  if  J5  ~  10”2  s“l.  In  such  a  case,  the 
recycling  of  Na02  to  Na  would  occur  via  reaction  RR5,  and  the  constraints  de¬ 
rived  here  for  reaction  RR6  can  be  easily  translated  into  similar  constraints 
for  reaction  RR5.  Photolysis  of  NaC>2  is  also  important  in  the  stratosphere, 
where  the  concentration  of  atomic  oxygen  is  small. 

The  value  adopted  for  reaction  RR4  by  Rowland  and  Makide  (1982)  assumes 
that  all  the  observed  EUV  absorption  of  NaOH  is  due  to  photolysis.  If  this  is 
not  the  case,  and  J4  is  smaller,  our  models  would  produce  more  NaOH,  its  day- 
to-night  change  would  be  smaller  and  the  effects  of  coupling  of  Na-NaO  to  the 
NaOH  reservoir  would  be  consequently  reduced. 

Ionization  of  atomic  sodium  can  also  affect  our  results  in  the  region 
above  the  peak.  Atomic  sodium  is  ionized  both  by  photons  and  by  change  ex- 
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Results  of  Models  A,  B,  C,  and  0  are  summarized  in  Table  4.  The  experi¬ 
mental  results  from  .Table  2  are  included  for  comparison.  Also  shown  in  Table 
4  is  the  midnight  D-line  emission  rate.  Results  on  the  D-line  emission  rate 
are  reasonable  for  all  models,  and  thus  are  inconclusive  for  our  purposes.  As 
noted  previously.  Models  C  and  D  agree  poorly  with  observations  of  the  peak 
height,  peak  density,  the  column  density,  and  its  diurnal  variation.  We  can 
therefore  rule  out  Models  C  and  D  as  unrealistic  descriptions  of  sodium  che¬ 
mistry  in  the  upper  atomsphere.  The  rate  for  RR6  (Na(>2  +  0  NaO  +  O2)  must 
therefore  be  significantly  lower  than  10_1°  cm3  s_1 ,  closer  to  10~13  cm3  s'1. 

It  is  harder  to  draw  any  definite  conclusions  to  differentiate  between 
Models  A  and  B.  The  best  basis  for  differentiation  is  the  diurnal  variation 
of  local  densities  and  total  column  densities  of  Na.  Local  density  variations 
for  both  models  are  consistent  with  observations,  given  the  altitude  resolu¬ 
tion  of  lidar  measurements.  Variations  in  the  total  column  density  for  Models 
A  and  B  are  more  sharply  differentiated,  and  the  question  thus  lies  on  the 
interpretation  of  the  experimental  data. 

Early  lidar  measurements  of  the  sodium  layer  yielded  day-to-night  ratios 
of  at  most  10%  (Gibson  and  Sandford,  1972)  or  changes  during  the  night  of 
about  4%  (Simonich  et  al.,  1979).  More  recent  observations  of  Clemesha  et  al. 
(1982)  which  were  carried  out  on  a  continuous  basis  yielded  diurnal  variations 
with  an  amplitude  of  ~  15%  from  the  mean  (30%  peak-to-peak).  The  observations, 
however,  exhibited  a  strong  semidiurnal  component,  and  the  above  authors  sug¬ 
gested  that  the  diurnal  variation  is  due  entirely  to  tidal  effects.  Thus, 
definite  conclusions  could  be  made  only  after  incorporation  of  tidal  effects 
into  existing  photochemical  models,  coupled  with  higher  resolution  measure¬ 
ments  of  the  diurnal  variation  of  Na.  We  thus  conclude  that  a  rate  of 
10  10  cm3  s  1  for  RR7  (NaO  +  H2O  +  NaOH  +  OH)  is  still  consistent  with  present 
experimental  data. 

Another  consequence  of  the  high  rate  for  reaction  RR3  is  that  the  diurnal 
variation  of  atomic  sodium  below  the  peak  can  be  estimated  to  first  order  in 
Models  A  and  B  from  the  photochemical  equilibrium  relation  between  Na  and  Na02 
given  in  equation  (5).  The  diurnal  variation  of  Na  then  follows  that  of 
atomic  oxygen  if  Na02  does  not  change  appreciably,  a  condition  satisfied  below 
the  Na  peak  in  Models  A  and  B.  Comparison  cf  the  diurnal  variation  of  0  in 


ratio  can  be  estimated  from  the  approximate  chemical  equilibrium  conditions 
for  Na02! 

k3[Na][02][M]  ~  k6[Na02][0].  (5) 

In  both  Models  A  and  B,  Na  is  the  most  abundant  species  at  90  km.  At  lower 
altitudes,  with  larger  values  for  [O2]  and  [M]  and  smaller  values  for  [0], 

Na02  becomes  the  most  abundant  species  at  82  km. 

In  Model  A,  with  production  of  NaOH  from  the  reactions  of  Na02  with  H  and 
OH  alone,  the  [Na0H]/[Na02]  ratio  is  approximately  given  by 

(Ji*  +  k12  [H] )  [NaOH]  =  [Na02]  (k10[H]  +  kn[0H]).  (6) 

Rapid  changes  in  [NaOH]  are  expected  at  dawn  and  dusk  due  to  the  onset  and 
disappearance  of  photolysis.  However,  the  concentration  of  NaOH  in  Model  A  is 
sufficiently  small  at  90  km  that  its  variation  does  not  affect  the  diurnal 
behavior  of  Na  and  Na02  (Figure  7b).  There  is  no  calculated  diurnal  variation 
for  Na  and  Na02  at  90  km,  because  0  shows  no  diurnal  variation  at  this  alti¬ 
tude  (Figure  2).  At  82  km,  the  calculated  Na02  shows  a  higher  concentration 
in  the  daytime  (Figure  7a)  because  of  the  interconversion  with  NaOH  which  has 
comparable  concentration  during  the  night.  The  diurnal  variation  of  Na  fol¬ 
lows  closely  that  of  0  as  suggested  by  the  relation  in  equation  (5). 

In  Model  B,  equation  (6)  has  to  be  modified  to  Include  an  additional 
source  from  the  reaction  of  NaO  +  H2O.  As  a  result,  the  calculated  [NaOH]  is 
much  higher  than  in  Model  A.  Because  of  Its  larger  concentration,  NaOH  is 
more  effective  in  influencing  the  diurnal  behavior  of  Na  and  Na02.  For  in¬ 
stance,  a  slight  diurnal  variation  is  seen  in  both  Na  and  Na02  at  90  km  (Fi¬ 
gure  7d).  At  82  km,  however,  the  diurnal  behavior  of  Na  is  still  primarily 
controlled  by  the  behavior  of  0. 

The  primary  difference  between  Models  A  and  B  lies  in  the  concentration 
of  NaOH  which  has  not  been  observed.  Unfortunately,  NaOH  exerts  only  a  se¬ 
cond-order  effect  on  the  diurnal  variations  of  Na  making  it  difficult  to  dif¬ 
ferentiate  between  the  two  models  based  on  Na  observations. 
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Noon  altitude  profiles  for  Na  are  shown  in  Figure  4  for  the  different 
models.  Measurements  by  Simonich  et  al.  (1979)  for  midnight,  and  Granier  and 
Megie  (1982)  for  noon  are  also  shown  at  altitudes  above  85  km.  As  indicated 
by  the  data,  the  diurnal  variation  of  Na  is  negligible  between  85  and  95  km. 
Agreement  between  these  measurements  and  models  C  and  D  is  poor.  Better 
agreement  is  obtained  with  Models  A  and  B,  especially  in  the  region  of  the 
peak.  The  discrepancy  between  observation  and  the  calculated  results  above 
95  km  is  probably  due  to  the  effect  of  ion  chemistry.  We  will  discuss  this  in 
more  detail  in  the  next  section. 

The  diurnal  variation  of  the  total  Na  abundance  is  presented  in  Figure  5, 
again  for  all  four  models.  The  absolute  magnitude  of  the  abundance  is  in  good 
agreement  with  observed  values  for  Models  A,  B,  and  C,  while  Model  D  yields 
column  densities  higher  than  those  observed  for  average  conditions.  The 
diurnal  variation  of  the  total  abundance  is  consistent  with  observations  for 
Models  A  and  B  (5%  and  15%,  respectively),  but  is  too  large  for  Models  C  and  D 
(30%  and  45%,  respectively).  The  large  diurnal  variations  calculated  for 
Models  C  and  D  follow. from  the  results  that  much  of  the  sodium  is  found  below 
82  km  where  the  diurnal  deviations  of  Na  is  the  largest. 

The  results  presented  in  Figures  4  and  5  suggest  that  Models  C  and  D  do 
not  reproduce  the  observed  Na  altitude  profile  and  diurnal  variation  of  the 
sodium  abundance.  Models  A  and  B  yield  very  similar  results.  Further  discri¬ 
mination  between  these  two  models  would  have  to  be  based  on  the  diurnal  be¬ 
havior  of  the  sodium  species  at  specific  altitudes. 

Figures  6a  and  6b  illustrate  the  diurnal  variation  of  Na  at  82,  84,  and 
86  km,  for  Models  A  (Figure  6a)  and  B  (Figure  6b).  Although  the  diurnal  vari¬ 
ation  in  Model  B  is  slightly  larger  than  in  A,  both  models  yield  results  in 
satisfactory  agreement  with  experimental  data.  We  should  note  that  the  day- 
to-night  ratio  changes  dramatically  over  an  altitude  range  of  2  km  in  this 
region.  This  fact  should  be  kept  in  mind  in  comparing  theoretical  results  to 
lidar  observations,  since  the  latter  typically  have  an  altitude  resolution  of 
1  -  2  km  (Megie  and  Blamont,  1977;  Simonich  et  al.,  1979). 

The  interaction  between  the  different  sodium  reservoirs  is  illustrated  in 
Figures  7a  through  7d.  The  results  can  readily  be  explained  in  terms  of  the 
reaction  rates  assumed  for  each  model.  In  the  region  80-90  km,  the  [Na]/[Na02] 
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Altitude  profiles  of  atomic  sodium  at  noon  calculated  for  models  A, 

B,  C,  and  D.  The  measurements  of  Simonich  et  al.  (1979)  at  midnight 
(x)  and  of  Granler  and  Megie  (1982)  at  noon  (o)  are  also  included  for 
comparison. 
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Figure  5.  Diurnal  variation  of  the  total  column  abundance  of  atomic  sodium, 
calculated  for  models  A,  B,  C,  and  D. 


